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The interaction of a magnetic vortex with a rotating magnetic field causes the nucleation of a 
vortex-antivortex pair leading to a vortex polarity switching. The key point of this process is the 
creation of a dip, which can be interpreted as a nonlinear resonance in the system of certain magnon 
modes with nonlinear coupling. The usually observed single-dip structure is a particular case of 
a multidip structure. The dynamics of the structure with n dips is described as the dynamics of 
nonlinearly coupled modes with azimuthal numbers m = 0, ±n, ±2n. The multidip structure with 
arbitrary number of vortex-antivortex pairs can be obtained in vortex-state nanodisk using a space- 
and time- varying magnetic field. A scheme of a possible experimental setup for multidip structure 
generation is proposed. 

PACS numbers: 75.10.Hk, 75.40.Mg, 05.45.-a, 85.75.-d 



I. INTRODUCTION 

A magnetization curling occurs in magnetic particles of 
nanoscale due to the dipole-dipole interaction. In partic- 
ular, the vortex state is realized in a disk shaped particle, 
where the magnetization becomes circular lying in the 
disk plane in the main part of the sample, which possesses 
a flux-closure state. At the disk center there appears 
an out-of-plane magnetization structure (the vortex core, 
typically from 10 nmP to 23 nnP) due to the dominant 
role of the exchange interaction inside the coreP^ The 
vortex state of magnetic nanodots has drawn much at- 
tention because it could be used for high-density mag- 
netic storage and miniature sensors Apart from that, 
such nanodots are very attractive objects for experimen- 
tal investigation of the vortex dynamics on a nanoscale. 

An experimental discovery of a vortex core reversal 
process by excitation with short bursts of an alternating 
field 6 initiated a number of studies of the core switch- 
ing process. The mechanism of the vortex switching is 
of general nature; it is essentially the same in all sys- 
tems where the switching was observed.™^ There are 
two main stages of the switching process: (i) At the first 
stage of the process the vortex structure is excited by the 
pumping, leading to the creation of an out-of-plane dip 
with opposite sign nearby the vortex. The appea rance 
of such a deformation is confirmed experimentally!^^ 
(ii) At the second stage, when the dip amplitude reaches 
the maximum possible value, there appears a vortex- 
antivortex pair from the dip structure. The further dy- 
namics is accompanied by the annihilation of the original 
vortex with the new antivortex, leading effectively to the 
switching of the vortex polarity; the dynamics of this 
three body problem can be described ana 
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While the second stage is well understood, the physical 
picture of the first stage of the switching process, which 
is a dip creation, is still not clear. 

The aim of the current study is to develop a theory of 
the dip creation, which is a key moment in the vortex 



switching process. The dip always appears as a nonlin- 
ear regime of one of the magnon modes. In most studies 
the dip appears by exciting a low frequency gyromode, 
which corresponds to the azimuthal quantum number 
m = — 1; the frequency of this mode ujq lies in the sub 
GHz range. Excitation of the gyromode always leads to 
a macrospopic motion of the vortex as a whole; moreover 
the vortex has to reach some critical ve locity v CY i of about 
300 m/s in order to switch its polaritjP^H 

Recently we have reported about vortex core switch- 
ing under the action of a homogeneous rotating magnetic 
field B = B x +iB y = B exp(icj£)P The theory of the 
dip creation was constructed very recently in our previous 
paper 20 , where we found a dip as a nonlinear regime of 
the high frequency mode with m = 1; the switching pro- 
cess for that case is not accompanied by a vortex motion 
at all. In this paper we show that the dip can be excited 
for any azimuthal mode m by using a nonhomogeneous 
rotating magnetic field of the form (16). 



The paper is organized as follows. In Sec. [II] we for- 
mulate the model and describe the approach of a rotat- 
ing reference frame. Two kinds of numerical simulations 
are presente d in S ec. |III| micromagnetic OOMMF simu- 
lations (Sec. Ill A) and s pin-la ttice SLASI simulations in 
the rotating frame (Sec. IIIB). An analytical approach 
is presented in Sec. [TV) We discuss our results in Sec. [V] 
In Appendix [A] we prove the conservation law for the to- 
tal momentum J z under the action of the magnetostatic 
interaction. A possible experimental setup for multidip 
structure generation is discussed in Appendix [B] 



II. MODEL AND CONTINUUM DESCRIPTION 

The continuum dynamics of the the spin system can 
be described in terms of the magnetization unit vector 
m = M/Ms = (sin$cos(/?, sin$sin(/?, cos$), where $ 
and (p are functions of the coordinates and the time, and 
Ms is the saturation magnetization. In the subsequent 
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text only disk shaped samples are discussed. Therefore 
it is convenient to introduce dimensionless coordinates 
p = (x,y)/L within the disk plane and £ = z/h along the 
disk axis, where L is the disk radius and h is its thickness. 
The magnetization is assumed to be uniform along the z- 
axis and so the corresponding angular coordinates of the 
magnetization read 0(p) = $(Lp) and (j)(p) = ip(Lp). 
The energy functional of the system under consideration 
consists of three terms, 



E = 4ttM|V {s e 



(1) 



Here V is volume of the sample. The dimensionless en- 
ergy terms are the following: 
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is the exchange energy with £ = \/ A/AirMg being the ex- 
change length, A being the exchange constant. The mag- 
netostatic energy <f ms comes from the dipolar interaction, 
see Appendix|Aj and in the continuum limit it can be pre- 
sented as a sum of three terms: <f ms = ^ v ms + ^ s ms + <f e ms - 
Here 
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is the energy of the interactions of volume magnetostatic 
charges, where e = h/L is the disk aspect ratio and 

Iv dV = Jo 1 d C Jo 1 d PP Jo* d * with X) bein g P olar 
ordinates within the disk plane. 
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is the energy of the interactions of charges on the upper 
and bottom surfaces, where J § dS = dpp d%- 
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is the energy of the interactions of edge surface charges 
with the volume charges. Here E is the disk edge sur- 
face and J E dE = J Q d£ J n d%. Due to magnetization 
uniformity along the z-axis the distribution of upper and 
bottom surface charges is antisymmetrical and the distri- 
butions of volume and edge surface charges are uniform 
along the z-axis. As a result the interaction energy of 
upper and bottom surface charges with volume charges 
as well as with edge surface charges is equal to zero. The 
last term in describes an interaction with a nonho- 
mogeneous rotating magnetic field, see below. 



The evolution of magnetization can be described by 
the Landau-Lifshitz-Gilbert (LLG) equation: 
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Here and below the overdot indicates derivative with re- 
spect to the dimensionless time 



r = coot, coo = Att'jMs, 



(7) 



where 7 is the gyromagnetic ratio, 77 is the Gilbert damp- 
ing constant, and the factor tt appears due to the disk 
volume normalization. These equations can be derived 
from the following Lagrangian 
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d 2 p(l-cos<9)(/>-^ 



(8) 



and dissipation function 

/ d V(fl 2 +sin 2 00 2 ). (9) 

Let us start with the no-driving case S { — 0. The 
exchange interaction ([2| provides the conservation of the 
total magnetization M z along the cylindrical axis z and 
the conservation of the z-component of the orbital mo- 
mentum L z 



M z 



i J d 2 p cos 0, L z = - i J d 2 p cos 6d x (j) (10) 



due to the invariance under rotation about the z-axis in 
spin-space and physical space, respectively. It is well- 
known 21 that the magnetostatic interaction breaks both 
symmetries. Nevertheless for thin cylindrical samples, 
when the magnetization distribution does not depend on 
the thickness coordinate, the magnetostatic energy is in- 
variant under two simultaneous rotations 

X^X + ^o, (11) 
leading to the conservation of the total momentum 

Jz = lj d 2 p (cos - 1) (1 - d x (j)) , (12) 

see the proof in Appendix [A} 

A. Planar vortex and magnon modes 

The ground state of a small size nanodisk is uniform; 
it depends on the particle aspect ratio e: thin nanodisks 
are magnetized in the plane (when e < s c ~ 1.812) and 
thick ones along the axis (when e > e c ). When the par- 
ticle size exceeds some critical value, the magnetization 
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curling becomes energetically preferable due to the com- 
petition between the exchange and dipolar interaction. 
For a disk shape particle there appears the vortex state. 

Let us consider the static solutions of the Landau- 
Lifshitz equation with in-plane magnetization (0 = 7r/2). 
The in-plane magnetization angle (j) satisfies the Laplace 
equation V 2 (/> = 0. Typically for the Heisenberg magnets 
the boundary is free, which corresponds to the Neumann 
boundary conditions. However, the dipolar interaction 
orients the magnetization tangentially to the boundary 
in order to decrease the surface charges. Effectively this 
can be described by fixed (Dirichlet) boundary condi- 
tions. The simplest topologically nontrivial solution of 
this boundary value problem is the planar vortex, situ- 
ated at the disk center 
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(13) 



where (£ = ±1 is the vortex chirality. Such planar vor- 
tices are known for the Heisenberg magnets with strong 
enough easy-plane anisotropy^. Qualitatively, when the 
typical magnetic length is smaller than the lattice con- 
stant, the pure planar vortex can be realized. In the 
magnetically soft nanomagnets like permalloy, there ex- 
ist out-of-plane vortices. The typical out-of-plane vor- 
tex structure has a bell-shaped form with a core size 
about the exchange length. 3 

Magnons on a vortex background can be described us- 
ing the partial wave expansion: 



COS = COS<9 V 0) + ^«m,n(T)/| m |, n (p)e 2mX , 



(14) 



where a m , n a* m>n and /3 m>n /?* m>n . In the 
linear regime a m>n (r) = exp(ii? m>n T) and /8 m>n (r) = 
exp(zi? m5n r). We consider here the planar vortex, where 
the spectrum of eigenmodes is degenerate with respect 
to the sign of m. Typical eigenfrequencies are presented 
in Fig. [I] The functions / and g obey the following nor- 
malization rule 
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(fm,n9m',n') = / fm,n(p)dm f ,n f (p)P^-P — ^m,m , ^n,n f 



(15) 

Below we assume that for each azimuthal number m it is 
sufficient to take into account only one radial wave with 
a certain radial index n. Therefore the summation over 
n will be omitted. 



B. Interaction with a field 

Let us consider the effects of a magnetic field. The role 
of the field is to excite spin waves on a vortex background. 
It is well known that the low frequency gyroscopical mode 




FIG. 1: (Color online) Eigenfrequencies of the magnon waves 
in a vortex-state Py disk (150 nm diameter, 20 nm thickness). 
The spectrum was obtained using micromagnetic modelling 
with a f ixed in-plane vortex core (for technical details see 
Sec. |III| ). The frequencies are normalized by ujq — AtvjMs 
(30.3 GHz), n denotes the radial wave number - the num- 
ber of roots of the functions f(p) and g(p). The separation 
of waves with different m and n was achieved using spatio- 
temporal Fourier transform technique. 



can be excited by applying a homogeneous ac field with a 
frequency ujq in the sub GHz range. The strong pumping 
of such a mode is known to cause the vortex polarity 
switching. 

The vortex switching phenomenon was studied re- 
cently by applying a homogeneous high frequency rotat- 
ing magnetic field. 20 Such a field excites the azimuthal 
mode with azimuthal number m = ±1; the frequency uj\ 
of such a mode lies in the range of 10 GHzP^ 

In the present work we study the vortex switching pro- 
cess by exciting higher magnon modes with higher az- 
imuthal numbers m. In order to excite such a mode we 
propose to consider the influence of a nonhomogeneous 
rotating magnetic field 



B = B x +iB y = B e 



i(/j,+l)x+iwt 



fj,ez. (16) 



Such a field distribution is chosen, because it directly 
pumps the magnon mode with m = p. Using Ansatz 
( 14 ) one can easily calculate the Zeeman energy 
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Here and below we use the normalized field intensity 
b = Bo/4ttMs and field frequency Q = uj/ujq. For the 
case /i < — 1 the field with structure (16) can be created 



experimentally, for details see Appendix p] 
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FIG. 2: (Color online) Vortex state dynamics under influence of the field (16) with \i — —3, .Bo^OmT and u;=6GHz. Disk 
radius L=150nm. Columns a), b) and c) correspond to the different moments of time after the field switching-on moment: 24, 
80 and 150ps, respectively. The top row of these columns illustrates distribution of the out-of-plane magnetization component 
(m z ). The in-plane component distribution is shown in the second row. Isolines m x — and m y — are shown as dashed and 
solid lines, respectively. In the bottom row the corresponding two-dimensional spatial Fourier transforms of the distribution 
m z (r,x) are presented. Each bar of a certain color corresponds to a Fourier amplitude F mn , where the azimuthal wave 
number m is specified by the horizontal axis and the radial number n is indicated on the bar. The subfig. d) demonstrates a 
vortex-antivortex pair in detail. The used magnetic field at t = is shown in the subfig. e). 



C. Rotating frame of reference 

In the laboratory frame of reference the total energy $ 
depends explicitly on time because of the Zeeman term 
(17), which contains the explicit time dependence. But 



using a transition into rotating reference frame (RRF) by 
the way of 

= + ^r, x = x + ^r, fj=-, (18) 

and using the invariance of the exchange and magneto- 
static energies with respect to the simultaneous rotations 
( 11 ), we obtain the time-independent energy in the RRF: 



S = <r 



(19) 



Here the Zeeman term <f f does not contain time explicitly 
and has the form 



(20) 



The dissipation function ([9| in the RRF reads 



d 2 p 0d x + sin 2 6(0^ - 1)0 



(21) 



III. NUMERICAL SIMULATIONS 

To investigate the magnetization dynamics of a vortex 
state nanodisk under influence of the external magnetic 
field (16), two kinds of simulations were used. 



A. Micromagnetic simulations 

The main part of the numerical results were obtained 
using the full scale OOMMF micromagnetic simulations 
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with the material parameters of Permalloy NigiFeig: ex- 
change constant A = 1.3 x 10 _11 J/m, saturation magne- 
tization Ms = 8.6 x 10 5 A/m, damping constant r] = 0.01. 
The on-site anisotropy was neglected and the mesh cell 
was chosen to be 2nm x 2nm x h. For all simulations the 
thickness h = 20nm. 

An example of magnetization dynamics induced by the 
field (16) with fi = — 3 is explored in Fig. [5] 

Initially we have a nanodisk in a vortex state which 
is a ground state. The vortex polarity was chosen to 
be negative. During the first tens of picoseconds of the 
field influence the generation of a magnon mode with az- 
imuthal number m = —3 is observed (see column a). The 
Fourier amplitude which corresponds to m = appears 
due to the out-of-plane core component of the initial vor- 
tex. Due to the pumping the mode with m = \i goes to 
the nonlinear regime: areas with sign(m 2 ) = sign(— fiuj) 
become localized which can be interpreted as dips for- 
mation (see column b). The number of dips is equal to 
For small field amplitudes such a multidip structure 
achieves some stationary regime and rotates around the 
vortex center forever. This phenomenon is described in 
detail below. Fig. [2] demonstrates the example of the 
magnetization dynamics for larger field amplitudes: the 
vortex-antivortex pairs are nucleated from the dips (see 
column c) and subfig. d). The subsequent dynamics of 
the vortex-antivortex pairs is rather complicated: the 
trajectory of a pair motion is not circular, pairs radiate 
magnons during the motion and self-annihilate. Then the 
process repeats periodically. The question about vortex- 
antivortex pair motion in the presence of an immobile 
central vortex is an open problem. But in this paper 
we focus our attention on the dips creation mechanism 
only. It should be also emphasized that the initial vortex 
in Fig. [2] is not pinned but it remains immobile during 
the dynamics. Using the simulations we found that such 
a vortex stability is observed only when the number of 
dips is greater than 2. The theoretical explanation of this 
phenomenon is an open problem. 



Using external fields of the form (16) with different 
fi one can obtain a multidip structure with an arbi- 
trary number of dips. This possibility is demonstrated 
in Fig. [3} 

It is important to emphasize the following properties 
of the multidip structure formation: (i) All dips have the 
same polarity which is determined by the sign of m only 
and the dip polarity does not depend neither on the po- 
larity of the central vortex (see the insets for \i = +3 and 
/i = +4) nor on the existence of an out-of-plane vortex 
core at all (see the insets for fi = — 2 and /i = —1). (ii) 
the multidip structure with n dips is composed by modes 
with azimuthal numbers m = n, 2n, 3n, . . . The contri- 
bution of a mode decreases when its azimuthal number 
m increases, (iii) Contribution of the mode m = is 
essential for a single-dip structure only (compare insets 
for /i = — 2 and // = —1). 

If the strength of the applied field is sufficiently small 
to prevent vortex-antivortex pair formation, the magneti- 
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FIG. 3: (Color online) Multidip structures created by the field 
(16) with different \i. To demonstrate the localized character 
of the dips the distribution of the out-of-plane component 
m z was obtained along the circle and the line £ which pass 
through the extremum points of a multidip structure. The 
distributions along the circle and the line £ are shown above 
Fourier spectra on the left and right plots, respectively. The 
other notation is the same as in Fig. [2] For the cases \i — —2 
and {i — — 1 the central vortex core is planar and artificially 
fixed, for the cases fi — +3 and \i — +4 the central vortex has 
opposite polarity and is not fixed. 



zation dynamics reaches some steady-state regime. This 
is demonstrated in Fig.|4]a) which is based on simulations 
with a fi = +3 field. Such a field creates potentially a 
structure with three dips of negative polarity. The time 
dependence of the minimal value of the m z component 
of the azimuthal mode with m = \i is shown. In case of 
multidip structures the plotted quantity is the depth of 
the dips, therefore the notation m^ ip will be used here- 
after. One can see that the value m^ ip reaches some 
steady-state level which depends on the frequency of the 
applied field. This dependence has an unusual resonance 
character, see Fig. [4]b). 

The transition from the regime of linear modes to the 
multidip regime occurs sharply when the frequency of the 
applied field reaches some value cjr es . Vertical steps on 
the dependencies m^ ip (uj) in Fig. [djb) correspond to the 
above indicated transition. It should be noted that the 
critical frequency depends on the field amplitude and is 
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FIG. 4: (Color online) Characteristics of the process of mul- 
tidip structure formation. The data was obtained from sim- 
ulations with a Py disk with radius L=76nm and thickness 
/i=20nm. If the multidip structure is formed, then m^ ip de- 
notes the depth of a dip, otherwise it denotes the amplitude 
of the corresponding azimuthal mode, a) - the process of 
steady-state three-dips structure formation for a certain field 
strength and for different frequencies, b) - resonant-type de- 
pendencies of the steady-state dips depth on the field fre- 
quency for different field amplitudes and for a certain fi. c) - 
the same as b) , but for a certain field amplitude and different 
\i. d) - data extracted from the resonant curves b): resonant 
frequency (left axis, solid line) and maximal dips depth (right 
axis, dashed line) vs applied field amplitude. 



always smaller than the eigenfrequency uj^ of the corre- 
sponding mode, see Fig. [4]d). This points to an inher- 
ently nonlinear nature of the resonance which is discussed 
in Section lTVl Similar resonances are a feature of the mul- 
tidip structures with different number of dips, Fig. [4]c). 



B. Spin— lattice simulations 

We have described above the micromagnetic study of 
the nonlinear dynamics of the magnetization under the 
influence of a time-dependent magnetic field. The main 
issue of this study is the creation of a multidip struc- 
ture, i.e. a stable nonlinear state of the system, which 
rotates due to the field rotation. Since the total energy 
of the system becomes time independent in the RRF, see 
(19), one can suppose that a multidip structure forms a 



stationary state of the system in the RRF. 

In order to check the RRF approach, we used another 
kind of simulations. Namely, we performed SLASI simula- 
tions, an in-house-developed spin-lattice code™ SLASI 
simulations are based on the numerical solution of the 
discrete version of the LLG equations ^ 



dS^ 

dt 



H 
s 



dt 



(22) 




FIG. 5: (Color online) Multidip structure formation under 
the influence of the field (16} with /i = -3 and 6=0.008. The 
data were obtained from SLASI simulations for a disk with 
radius L = 50a, thickness h = 5a, exchange length £ = a, and 
damping coefficient 77 = 0.1. a) Solid curves correspond to the 
simulations in the laboratory reference frame for the rotating 
magnetic field with different frequencies The dashed curve 
corresponds to the simulations in the RRF with Q rot — Q 
under the influence of a static field, b) Out-of-plane spin 
distributions from simulations in the RRF with Q rot — 0.3 at 
the moment r ~ 24 (same color code as in Fig. [2]). 



square lattices of size (2L) 2 ; the lattice is bounded by a 
circle of radius L on which the spins are free. The Hamil- 
tonian <K = lK ex + IK dip + !K f is given by the Heisenberg 



exchange Hamiltonian ( Al ), the dipolar energy (A3), and 
the field interaction energy IK f , which is the discrete ver- 
sion of ( pTrj ) . The 4th-order Runge-Kutta scheme with 
time step 0.01 /N z was used for the numerical integration 
of Eqs. (|22|). 



First of all we performed SLASI simulations in the lab- 
oratory frame of reference. Here the spin-lattice simula- 
tions agree with our micromagnetic results. In particular, 
the resonance behavior of the multidip structure is well- 
pronounced in Fig. [5|i) (solid curves), with a maximum 
dip amplitude about the frequency Q = 0.3. 

The total energy of the magnet in the continuum limit 
in the RRF (19) has two main differences in compar- 
ison with the laboratory reference frame. First of all, 



instead of the time dependent Zeeman energy (17), one 
has the influence of a constant field with the same inten- 
sity. Apart of this, there appears an additional rotation 
energy <o rot = QJ Z . For the discrete system one can per- 
form a similar transformation, using the discrete Zeeman 
energy of the interaction with a constant field. The dis- 
crete analogue of the rotation energy is determined by 
the discrete version of the J z momentum (12). Results 



where the 3D spin distribution is supposed to be inde- 
pendent of the z coordinate. We consider Eqs. (22) on 2D 



of numerical simulations in the rotating frame of refer- 
ence are plotted by the dashed curve in Fig. [5^t) for the 
rotational frequency Q rot = 0.3; they are in a good agree- 
ment with simulations in the laboratory reference frame 
for Q = 0.3. Note that we presented on Fig. [5^l) simula- 
tions in the RRF for small enough rotational frequencies. 
The reason is that J z is not a good characteristics of the 
discrete system: both exchange and dipolar interaction 
break the conservation of J 2 , since there is no rotation 
invariance of the discrete system. Therefore simulations 
in the RRF works well when the contribution of the J z 
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term in the discrete Hamiltonian is small enough, and 
the discreteness effects are small. This is valid for small 
enough frequencies of the rotations, see Fig. |5|i), or for 
high frequencies but small enough times, see Fig. [5b). 



IV. ANALYTICAL DESCRIPTION OF THE DIP 
CREATION 

In order to describe the problem analytically, we con- 
sider the problem in the RRF. As we have seen using 
SLASI simulations, the dip can be considered as a st a- 
tionary state of the system in the RRF. 

To gain some insight how the interaction with the mag- 
netic field (which is static in the rotating frame) together 
with the rotation provides the multidip creation we use 
the Ansatz ( fl4| ). In order to simplify the model we con- 
sider the dip formation on a background of a pure in- 
plane vortex with cos# v =0. This approximation is con- 
firmed by our numerical simulations, as well as by a pre- 
vious studyP^ The field pumps directly only the mode 
with azimuthal number m = ji. This mode is coupled, 
first of all, to the mode with m = — /i, 0, ±2/i. Thus we 
will consider modes with m = 0, ±2/i. Taking into 
account only cubic nonlinear terms in the magnetic en- 
ergy, one can use the Ansatz (14) to calculate the effective 
Lagrangian 



m=0;±/x;±2/x 

and the effective dissipation function 

2imQ 



(23) 



j? eff = 1 



E A 



/i 



+'B m /3 m ( f3 m f3 m 



(24) 



where A m (f^) and £ m = (g^). 

The effective energy consists of several parts 
+ Cff + + where 



(25) 



describes the linear part of the modes oscillation, the 
energy of the interaction with the field <^ e f ff is equal to S { 
in (20), the rotation energy 



(26) 



appears due to the transition into the noninertial frame 
of reference, and the energy of the nonlinear coupling 
between the modes has the form 



ra+n ' 



(27) 



For the explicit form of the magnon frequencies i? m and 
the nonlinearity coefficients k^ ■ see Appendix [cj 

Assuming that the nonlinear coefficients take approxi- 
mately the same values fc| . w k one can obtain the equa- 
tions of motion for the amplitudes a m , /3 m in the form 



m 



/i 



-ik ^ m'a^OLm'+m + T]?> m I $ m + i — flfim ) , 
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(28) 



In the infinite set of equations (|28| we restrict ourselves to 
the finite number of equations which do not contain any 
amplitudes a^ft, with the exception of amplitudes with 
indices i = 0, d=/x, ±2/i. The obtained system was solved 
numerically for the case \i = 3 and different field ampli- 
tudes. The values of the corresponding eigenfrequencies 
^o,/i,2/i were chosen to be equal to the ones obtained from 
extra micromagnetic simulations of the magnon dynam- 
ics, these values are marked by dashed circles in Fig. [I] 
The frequencies Qq and Q\^\ correspond to the lowest 
modes with radial number n = while corresponds 
to the mode with n = 1, because of two reasons: (i) ac- 
cording to the Fourier spectra in Fig. [2] and Fig. [3] the 
mode with this radial number dominates among modes 
with azimuthal number m = 2/i; (ii) as it will be shown 
later the second resonance R 2 (see Fig. |I}b) appears at 
the frequency i?2u/2 which corresponds to the mode with 
n = 1 (see Figjl]). According to the mentioned spectra 
the amplitude of the dips in a multidip structure is deter- 
mined mainly by the value a M , i.e. the amplitude of the 
out-of-plane component of the mode with m = fi. The 
numerically obtained dependencies are shown in 

the Fig. . These dependencies are in good agreement 
with the behavior of the dips depth, obtained using simu- 
lations, see Fig. |4^i). Moreover, the steady-state value of 
a M has the same properties as the steady-state dips depth 
m^ ip : (i) its frequency dependence has a resonance char- 
acter, (ii) the resonance frequency is lower than the eigen- 
frequency, see Fig. [6^i). Amplitude-frequency character- 
istics obtained numerically from model (28) (Fig. |6]>d) 
let us conclude that the phenomenon of an abrupt ap- 
pearance of a dip (multidip) structure in a vortex state 
is a nonlinear resonance in a system of nonlinearly cou- 
pled modes with m = 0, ±/i, ±2/i. The corresponding 
resonance transition is shown as R\ in Fig. |4Jd) and in 
Fig.|6b)-d). 



In the no-damping and weakly nonlinear limit (k <C 
1) the stationary values a m and /3 m can be represented 
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FIG. 6: (Color online) The nonlinear resonance in system 
(28). Subfigure a) demonstrates time dependencies of the 
amplitude of the out-of-plane component of the mode with 
m = fjb. The shown dependencies were obtained numerically 
directly from ( |28| ) for different field frequencies Q (the values 
are in the oval frames). The values of the other parameters 
were the following: ji = 3, = 0.4, J? = 0.47, J? 2/ x = 0.92, 
b 0.01, J] 0.02, k 0.1, Ai = ¥>i = 1. Amplitude- 
frequency characteristics for different modes are presented in 
the subflgures b)-d). The numerical values of the parameters 
are the same as in subfigure a) . The thick (red) curve is build 
for the field amplitude b = 0.01, thus it fully corresponds to 
the subfig. a). The small insets shows the mentioned curve 
in full scale. 



can be strictly described as the dynamics of nonlinearly 
coupled modes with azimuthal numbers m = 0, ±n, ±2n. 
The multidip structure with an arbitrary number of dips 
(or vortex-antivortex pairs) can be obtained in a vortex- 
state nanodisk using a space- and time- varying magnetic 
field of the form ( [l6| ). A scheme of a possible experimen- 
tal setup for multidip structure generation is proposed in 
Appendix |B) 
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Appendix A: Dipolar interaction and the 
conservation of the total momentum 



approximately as follows: 
ib e Q 



Q 2 ' 



a 



n 2 

2kfiQ^ \a^\ 



-2knb 2 e f2^n 



fc/ua 2 2f2 2 + Q^Qip 



.2 ' 



O 4ft 2 



fan = ikfMXp 



2fl„ 



2fi 



4ft 2 



(29) 



According to the last two equations an additional res- 
onance is expected for the frequency Q = i?2^/2. This 
resonance is observed in the simulation results (see Fig.|4j 
b) and it is well recognized in Fig. [6] b-e, is denoted as 
R 2 - 



V. CONCLUSIONS 

We have presented a detailed study of the dip struc- 
ture generation, which always precedes the vortex polar- 
ity switching phenomenon. The physical reason for the 
dip creation is softening of a magnon mode and conse- 
quently a nonlinear resonance in the system of certain 
magnon modes with nonlinear coupling. The usually ob- 
served single-dip structure is a particular case of a multi- 
dip structure. The dynamics of the structure with n dips 



In this Appendix we consider a ferromagnetic system 
described by the classical Heisenberg isotropic exchange 
Hamiltonian 



(Al) 



and the dipolar interaction !K dip : 



:K dip 



D 

~2 



£ 

n.n' 



•3(5,, 



>)(S n 



,/|3 



(A2) 

Here S n = (5^, 5^, S^) is a classical spin vector with 
fixed length S in units of action on the site n = 
(n x , n y: n z ) of a three-dimensional cubic lattice with inte- 
gers n x , n y , n z , J is the exchange integral, the parameter 
D = 7 2 /a 3 is the strength of the long-range dipolar inter- 
action, 7 = g\e\/(2mc) is the gyromagnetic ratio, g is the 
Lande-factor, a is the lattice constant; the vector 5 con- 
nects nearest neighbors, and e nn > = (n — n') / \n — n'\ 
is a unit vector. 

Our main approximation is that S n depends only on 
the x and y coordinates. Such a plane-parallel spin dis- 
tribution is adequate for thin films with a constant thick- 
ness h = N z a and nanoparticles with small aspect ratio. 
Using the above mentioned approximation the dipolar 
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Hamiltonian can be written as follows 
IK dip = — — ^2 \Avv> (S u • S U f — 3S*S*,) 



(A3) 



Here the sum runs only over the 2D lattice. All the in- 
formation about the original 3D structure of our system 
is in the coefficients A^/, B uu f and C^/, 



2 ^ 

n z ,n z 

3 ^-^v 

2 f-^ 

n' z ,n z 

3 E : 

n z ,n z 



a 2 

ynn' 



(A4) 



rirv ynn' 



where we used the notations: x Tl 



a(n y 



= a(n x - n' x ), 



a(n z - n' z ), p„„> 



\/ X nn' + Vnn f ' r ™™' — \/ ptv' + ^n/ " 

The continuum description of the system is based on 
smoothing the lattice model, using the normalized mag- 
netization m(r) = (g/^B /a 3 Ms) J2 n S n 5(r — r n ). Then 
the exchange energy <f ex (normalized by 47rMj), the con- 
tinuum version of ( |A1| ), takes the form The normal- 
ized magnetostatic energy, which is the continuum ver- 
sion of (A3), is 



1 

16tt 



d 3 r/dV W[m ' m1 



W[ra,ra'] = (R 2 -3| 



— 2 cos 6 cos (9' 



z — z 



sin (9 sin 0' cos (</> — <//) 
p 2 cos (0 + 0' - 2x) 



3 sin # sin 0' 

+ // 2 cos (0 + </>' - 2x0 - W cos (0 + 0' - X ~ X f ) 
R(r,r') - 



p 2 + p' 2 



2pp' cos(x - xO + (z ~ z') 2 - 

(A5) 



The magnetostatic energy in the form (A5) is invariant 



under to simultaneous rotations of <p and x with the same 
constant angle ipo, see Eq. (11). 



The consequence of such an invariance is the conserva- 
tion of the total momentum (12). Let us show explicitly 



that J z is conserved. The time derivative of the total 
momentum 



dJ z 
dt 



d 6 rd y 



l-cos(9 dS 
sin ' ~d6 



I 



d 3 r(cos0- 1) [d x ,d t ] <j>, 



(A6) 



3 # SS 53 
d r ( — - ^r-o Y 



where we used an explicit form of Eqs. Q in the case 
of absence of magnetic field and damping. The first 
term in ( A6 ) vanishes due to the cylindrical symmetry of 
the sample. The second term contains the commutator 
[d X: d t ] </>; it can take nonvanishing values for the singu- 
lar field distributions like 2D solitons and vortices with 
— QX + const, which results in [d x ,d y ]4> = 2irq5(r)^ 
Neverthe less this singularity does not influence the last 
term in ( |A6[ ) due to the vanishing factor (cos0 — 1) at 
the singularity point. This replacement cos0 — > cos0 — 1 
corresponds to the regularization of the LagrangiarP3. 
Let us discuss the last term in (A6). Since an isotropic 



exchange interaction allows the conservation of M z and 
L z separately, we need to discuss here the influence of the 
mag netostatic interact ion only. Using the explicit form 
(A5), one can rewrite (A6) as follows 



dJz_ 
dt 



1 fdVdV 
8tt/ R 5 



dW dW 



h + h 

8tt 



Here the contribution 1\ takes the form: 



d / W 



-W 



d 



d x \R 5 



The first integral vanishes due to periodicity on x an d 
cylindrical symmetry. The derivative in the last term 

%UJ = "^ rSm(x " x) 

is asymmetric with respect to the replacement r r', 
hence this term vanishes after the integration, and I\ = 0. 

The integral h = J d 3 rd 3 r'F[r, r'\ has the asymmet- 
rical kernel 

F[r,r'] =sin(9sin(9 / [(^ 2 -3|z-z / | 2 )sin(^-^ / ) 



p' 2 sin ( 



2x') 



p 2 sin ( 



-F[r',r], 



therefore, I 2 =0. 

Finally, one can state that the total momentum J z 
is conserved for a cylindrical sample under the action 
of magnetostatic interaction. One should note that the 
conservation of J z is known for the local model of the 
magnetostatic interaction <f™ c s = x J d 3 r (V • m) 2 P2 



Appendix B: Experimental possibility of a multidip 
structure creation 

Let us consider a set of long conductive wires orien- 
tated perpendicular to the disk plane. Let all these wires 
be uniformly spaced along a circle which is concentric 
with the disk and has radius a R > L. The described 
configuration is shown in Fig. 7] a). 

Let the n-th wire conduct the current I n = I cos (jut + 
27T77m/iV), where N is the total number of the wires, m is 
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FIG. 7: (Color online) Possible implementation of an exper- 
imental facility for pumping of a selected azimuthal magnon 
mode. Subfig. a) demonstrates t he u tility assembling. Sub- 
fig, b) shows the magnetic field (Bl) within the disk area. 
The number of the wires is N = 12 and the phase shift of the 
currents is = ty/2 (corresponds to m = 3, see text). Blue 
and red lines in subfig . c) illustrate dependencies m^ ip {t) ob- 
taine d us ing the field (16) with ji = — 3 (see Fig.[2]e) and the 
field (B2) with m — 3, respectively. In both cases the field 
amplitude was 15mT at the disk edge. L — 76nm, h =20nm. 



the azimuthal number of the mode we aim to excite and uo 
should be close to the corresponding eigenfrequency Q m . 
The described set of wires produces a magnetic field of 
the form 



n=0 



(r - R n ) 2 



(Bl) 



where R n = R[cos(27rn / N) , sin(27m/7V), 0] is the radius- 
vector of the n-th wire, (r, x, z) are coordinates of a cylin- 
drical frame of reference with e z being oriented directly 
to the reader, and c denotes the speed of ligh t. F or the 
case N ^> m one can proceed from the sum in (Bl ) to an 



integral. And for the case L <C R the obtained integral 
yields 



t\m-l\ 



NI r 
2Rc\ 

sm[ut - (m - l)x] [f h 
sm[ojt - (m + l)x] [^ |m+11 ■ 

^X) = "^{ 
cos[ujt-(m-l) X ] ^ |m_11 
cos[ojt — (m + l)x] 



^|m| 
. £|m| + l 



1 



t\m+l\ 



^|m| + l 
, £|m| + l 



}■ 



}■ 



(B2a) 



(B2b) 



where ( = r/i? « 1. The field (B2) corresponds to the 
field (16) with p = — |m| and sgn(u7) = sgn(m), but con- 



trary to (16) it is radially dependent. The radial depen- 



dence does not change the ability of the field (B2) for the 



production of multidip structures. The radial dependent 



field results only in a small decreasing of the dips ampli- 
tude, see Fig |7]c). The accuracy of conformity of (p2fo 
with (Bl) strongly depends on N. Using the minimal 



number of the wires N — 2m one can excite the standing 
wave with azimuthal number m. To obtain the multidip 
structure one should take 2m < N < 4m. It is senseless 
to take N > 4m because it does not improve the accu- 
racy significantly. Satisfactory accuracy was achieved for 
L/R < 1/2. 



Appendix C: Microscopic expressions for magnon 
frequencies and nonlinearity coefficients 



Substitution of the Ansatz ( 14 ) into Q and integration 
results in: 



■ = l^ 2 i;{\<*rn\ 2 (fZ + (rn 2 -l) fl 

m * 

\P m \ 2 (g' 2 



12/ -'2 ,. m 2«a 
p 2 



/o-^ 2 2 n-k I f rrif n Qm + n \ , n / 4 , \ . 
^ h? S Z^ a m a nPm+n m \ j>. ) + 0(a , /3 ). 



(CI) 



Using that for a weakly excited in-plane vortex state cp = 
X + £§ + ft, vn z = m z the divergence has the form 

P I dp dx V d Xj\ 

(C2) 

one can perform a similar action to obtain the corre- 
sponding magnetostatic terms: 



^s = ^J2\^\ 2 (7m(phm(p')) 



ms(v) 



IS 

71 



E 



m 



(p) 



Jn(p')f m (p')\ mS(V) 



(p) 



K g n {p')g m {p') 
p' 



P I m+n 

1 \ \ ms(v)" 



m+n 



^s = Y £ J2\ a m\ 2 (fm(p)fm(p')) 



ms(s) 
m ' 



+ 0(a 4 ,/? 4 ), 
(C3) 

(C4) 



IS 

71 



E 



m 



„ „ o* I fn(p)fn(p)9m+n(l) \ mS{e) 
QLn&rnPrn+n \ p J 



m+n 



iooo* 1 9n(p)gm(p)gm+n(l)\ mS{e) 
1 PnPmPm+n 



m+n 



+ 0(a 4 ,/? 4 ). 
(C5) 
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The notation 



, , _ dg m (p) g m (p) 
ImW- dp + z p 



is used here and three different types of averaging are 
defined: 



(C6) 



(F(p,p>))™ s{v) = (~l) k J dpp J dp'p'F(p,p>) J e ~ EX £2 ^ +£X J k (px)J k (p'x)da 



1 1 oo 

(F(p,p'))™ s(s) = {-l) k J dpp J dp'p'F{p,p') J {I- e-^)J k {px)J k {p'x)dx 







{F{p))T ie) = (~ir i dppF(p) 



1 + ex 



J k (px)J k (x)dx. 



Comparing (C3)-(C5) with (25) and (27) one can conclude that the eigenfrequencies read 
ft» =te 2 ( C + (m 2 - + - (fm(p)fm(p'))Z S(s) 



h 2 \ "" ' p* I e 

f 2 I a 2 

~h 2 



ms(e) 
m ' 



and the nonlinearity coefficients have form 



~l,oi _ A /o^ 2 Jl I fmfnQm+n \ g 



lm+n(p)- 



fn( P ')f m (p') X" 1 ^ / fn(p)fn(p)9m + n(l) 



m+n 



m+n 



fc' 3 =— 



lm+n{p)- 



i(p')g m (p') 

P' 



I \ \ ms(v) 



m+n 



gn(p)gn(p)gm+n(i)\ mS{e) 



m+n 



r 



(C7) 



(C8) 



(C9) 



It is interesting to emphasize the following properties of magnetostatic energy appear due to the volume charges 
the energy expansions: (i) The linear part of the mag- only; (iii) The exchange produces nonlinear terms of the 
netostatic energy proportional to \f3\ 2 is produced by the form aaf3. The same is true for magnetostatics, but ad- 
volume charges, while the part proportional to \a\ 2 is pro- ditional terms in the form /S/3/3 also appear, 
duced by surface charges; (ii) The nonlinear terms of the 
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